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ABSTRACT 


In view of the complex flight operation of the Space Shuttle propul- 
sion system together with an expected launch rate increase, the flight 
performance reconstruction process needs to be performed by automated 
computer programs. These programs must have the capability to quickly 
and reliably determine the true behavior of the various components of 
the propulsion system. For the flight reconstruction, measured values 
from the solid rocket motors, liquid engines, and trajectory are appraised 
through the Kalman filter technique to identify the most likely flight 
propulsion performance. 

A more detailed data collection program for the single SSME engine 
captive test firing evaluation is scheduled for startup in September of 
1988. Engine performance evaluation for the captive test firing requires 
a reconstruction process that is similar to the process that is used for 
the flight reconstruction. This paper describes analytical tools that may 
be used to reconstruct a propulsion system's true performance under 
flight and/or test conditions. 
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1 . 


INTRODUCTION 


During the summer of 1987, Rogers Engineering and Associates (REA) 
submitted a technical report to Marshall Space Flight Center. This report 
summarized the results of the Propulsion Estimation Development Verifica- 
tion performed by REA under contract NAS8- 36152. During that contract 
period, REA modified an existing program, which was developed under 
a previous contract , to include improved models of the Solid Rocket 
Booster, the Space Shuttle Main Engine (SSME) gain coefficient model, 
the vehicle trajectory using quaternions, and an improved Kalman filter 
based on the U-D factorized algorithm. In this report, this modified 
program will be identified by the symbols PROG1. 

Under a current contract, REA has proposed to design, evaluate, 
and refine a model for a single SSME that does not include the influences 
of other Space Shuttle components. The refined SSME model is expected 
to permit the collection of high quality measurements and to provide an 
improved quality of SSME performance estimation. REA's single SSME 
model should be available in the Fall of 1988, and this report refers to 
that program by the symbols PROG 2. 

Propulsion performance estimation procedures applied in programs 
PROG1 and PROG 2 are basically the same. Each requires the assumption 
that a correct model of the Space Shuttle's propulsion system exists and 
the model involves two main equations. They are an equation that models the 
dynamics of a state vector and a second equation that defines the rela- 
tionship between an observation vector and the state vector. Also, the 
engine performance estimation consists of using collected data on the 
observation vector to obtain a corresponding estimate of the state vector. 
This state vector estimate is called an estimate or a reconstruction of the 
engine's performance. 

The state vector estimate is based on the satisfaction of a statistical 
optimization criterion. Furthermore, data quality determination and the 
quality of the estimated state vector are judged through an application 
of statistics tools. Hence, in order for the Propulsion Laboratory to 
properly apply and correctly interpret the results of programs PROG1 
and PROG 2, an intuitive review of estimation procedures utilized by these 
programs is needed. 

Therefore, the objectives of this paper are to: 

1. Examine the assumptions and limitations of the Kalman filtering 
technique for the intended application to flight performance reconstruc- 
tion of propulsion systems and single engine performance reconstruction 
based on static testing data. 

2. Select specific topics from the area of statistics and probability 
for discussion and better understanding with EP55 personnel. 
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3. Identify some additional applications of contractor’s programs. 

Throughout this paper, underlined capital letters are used to denote 
vectors, capital letters denote matrices, and the identity matrix is denoted 
by the capital letter I. All vectors are of the column type and the 
transpose of any matrix or vector is denoted by using the letter T at 
the superscript position. The letter E denotes the expectation operator, 
N(U,Z) denotes the multivariate normal probability distribution with mean 
vector IJ and covariance matrix Z. The caret symbol written directly 
above a scalar or vector denotes a statistical estimator of that scalar or 
vector . 


2. NOTATION 


For any Space Shuttle flight or single SSME test, let denote the 

observed values of vector at time t. Each component of represents 

a relevant measurable output of the Space Shuttle Propulsion System or 
a relevant navigation measurement. For example, the components of 

may be oxygen pressurant flow, fuel volume flow, hydrogen pressurant 
flow, fuel flow pressure, fuel flow temperature, and chamber pressure 
whenever Z^ represents the output of a single SSME static test. For an 

actual flight of the Space Shuttle, Z t will contain all of the aforemen- 
tioned measurements for each of the three main engines plus additional 
components that represent navigational measurements and one measure- 
ment for each of the two Solid Rocket Motors (SRM). For the actual 
flight, Rogers (1987) lists 35 components for vector Z^ and 71 compo- 
nents for state vector X^. where is defined in the next paragraph. 
However, for the single engine static test setup, vectors Z^ and X t have 

fewer components because the SRM’s and navigational components become 
inactive . 

Vector X^ is a state vector of parameters to be estimated at time t. 
It is assumed that the observation vector is a function of the state 
vector X t> That is, 


Zj = h (X t ,t) + V t 


( 2 . 1 ) 


where h is some function and V t v N(0, R^). State vector X t is known 
to change with respect to time according to the equation 


x_t = f (x t ,t) + w t 


( 2 . 2 ) 
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where f is some function, W t v N(0,Q t ) and E[W t V t T ] = 0. 


If it is assumed that f and h in equations (2.1) and (2.2) are 
linear, then numerical procedures allow us to transform these equations 
into the form 


-k $ (k,k-l) ^k-1 + ^ 


(2.3) 


*k = H k X k + Yk 


(2.4) 


where k represents discrete values of t, ^ N(0,Q k ), v N(0,E k ) 
and * 

E(W k Vj. T ) = E(W i W. T ) = E(V i Vj T ) = 0 , 


provided i # j. The symbol $(k,k-l) denotes the transition matrix that 
connects vector to X ^ _ 

After a Space Shuttle flight has taken place or after a static test, 
the post flight /test reconstruction procedure seeks to use the observed 
values of and equations (2.3) and (2.4) to reconstruct (estimate) the 

vector X^ so that the estimation error is minimized. The reader should 

see Gelb (1974) for a review of minimization procedures. The estimated 

A 

value of X^ is denoted by X ^ . 


3. THE KALMAN FILTER SETUP 

Equations (2.3) and (2.4) along with assumed properties of vectors 
W j and allow for the development of the below Kalman filtering 

process. During the process, two types of estimators of X^ are possible 

for each k. A notation used by Gelb (1974) allows the two estimators 
to be distinguished. That notation is 

= the estimate of X^ using all observations up to and 
including observation (k - 1) . 

= the estimate of the error covariance matrix for X^ using 
all observations up to and including observation (k-1). 


X (_) 
2k 


(-) 
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X k = the estimate of using all observations up to and includ- 
ing observation k. 

P ^ = the estimate of the error covariance matrix for X, using 
k — K 

all observations up to and including observation k. 

The values of X, v , and P. v ' are obtained by using equation (2.3). 

The X. v ’ vector is often called the projected ahead value of X, and 

A (+) K 

P k v J is the projected ahead variance. Vector X k ' is called the up- 
dated estimate of X^ and P^ denotes the updated variance. Computation 
equations (3.1) through (3.5) describe the discrete Kalman filtering 


process. 


X k ( ‘ ) = *(k,k-l) x£i 

(3.1) 

P k (_) = *(k,k-l) P^ <D T (k,k-l) + Q k l 

(3.2) 

K k = p k ( ° H k T <H k p k ( ° H k - V 1 

(3.3) 

!k <+) = + K k <Zk - H k ik <_) > 

(3.4) 

V = « - K k H k> p k (_) • 

(3.5) 

A (-) 

To start the computation procedure, initial estimates 5 £q 

and P Q ( ^ 


are needed and are usually determined by the users of equations (3.1) 
through (3.5). 


The linear system described in equations (2.3) and (2.4) is essen- 
tial to the development of the discrete Kalman filtering process given in 
equations (3.1) through (3.5). Therefore, when functions f and h are 
nonlinear, as they are in the case of the flight /test reconstruction model, 
vector X^ has to be estimated through a linearization process. Two 

different linearization procedures are described in the next paragraph. 
The first procedure is called the linearized filter and the other is known 
as the extended Kalman filter. 
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Let be some known vector such that X^* = f (X^*,k). Often 

times X^* is called a reference solution or a reference nominal trajectory. 

That is, X^* is a known solution to equation (2.2) where the influence 

of Wj is not considered. If a discrete solution to equation (2.2) is X^ 

where X^ = X^ * + AX, then the linearization process may be used to 

estimate AX. Here, we see that for any AX and for any k value of t, 
a first degree Taylor series approximation of functions (2.1) and (2.2) 
may be determined by 



+ AX % f (X k *,k) + 



AX + W k 


and 


"3h-| 
- 9XJ 


CX k >k) + |_j . AX + ^ . 

— -k 


Since X^ = f (X k ,k), equations (3.6) and (3.7) reduce to 



(3.6) 


(3.7) 


(3.8) 


(3.9) 


Equation (3.8) is called the linearized dynamics equation and (3.9) is 
the linearized measurement equation. These two equations are linear and 
are equivalent to the linear equations (2.3) and (2.4). Hence, for each 
discrete time k, AX^ can be estimated and error covariance matrices can 

be determined. The state vector estimate at time k is then given by 



+ AX k 


( + ) 


(3.10) 
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where AX k is computed by applying the Kalman filtering process to 

equations (3.8) and (3.9). When each X k+1 (k = 0,1,..., T) is known 

prior to the beginning of the application of the Kalman filtering process 
to equations (3.8) and (3.9), the estimate in equation (3.10) is called 
the linearized Kalman estimate. When each X k+1 is determined from the 

previous estimate of X^, i.e., X * + ^ = & k + , the estimate in equation 

(3.10) is called the extended Kalman filter estimate. 


4. FAILURE DETECTION 


The current repalcement cost of a single SSME is about $50 million. 
This cost, combined with analyses of data from static test failures and 
shutdowns, suggest that there is a need for a more advanced state test 
failure detection system for the SSME. A 1985 report by Taniguchi 
emphasizes the importance of SSME failure detection improvements and 
identified some possible analytic designs for a failure detection system. 
Each of Taniguchi's failure detection designs assumed that a correct 
model of the single SSME exists. Each design also required the applica- 
tion of a statistics test as a decision tool. Therefore, a natural exten- 
sion of REA's SSME static test model would be into the area of failure 
detection during static testing. REA's static test model is already 
undergoing refinements and is expected to be available to Marshall Space 
Flight Center by the Fall of 1988. Also, according to Taniguchi (1985), 
data is available on 1200 static test firings of a single SSME. A brief 
overview of a failure detection setup is provided in the next paragraph. 

Assume that equations (2.3) and (2.4) are modified to yield 


X^ = 4>(k,k-l) + Aj (k,k-l) b k + W k 


(4.1) 


^k = 


H k*k 


+ A 2< k > b k + Vfc 


(5.2) 


where vectors bj^ are called biased vectors. If biases are regarded as 
additional states with dynamics such that 


^k+l 



then equations (4.1) and (4.2) may be transformed into the form 
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(4.3) 


Xk = 


$(k,k-l) A^k.k-l) 
0 I 


Xk-1 + 


I 

Lo. 




and 


2k = i H k A 2< k >i y k + Y*. 


where 


Xk = 


*k 

L^k J 


(4.4) 


If Kalman filtering is applied to equations (4.3) and (4.4), the bias 
vector b^ is estimated along with the other components of the state 
vector. 


If 

that bj^ 


no failure has occurred at time k, then it is intuitively reasoned 
= O. Of course, the judgement as to whether b. is really zero 


is determined by a multivariate statistical test. 


For additional details on bias vector estimation, the reader should 
consult Friedland (1983). 


5. CONCLUSION 


The assumptions and limitations of the Kalman filtering process have 
been discussed with NASA personnel. These discussions included rele- 
vant topics from the area of statistics. In particular, it was pointed out 
in this paper that equations (2.3) and (2.4) are the essential ingre- 
dients for the discrete variable reconstruction. Hence, if functions f 
and h of equations (2.1) and (2.2) are nonlinear, X^. can be estimated 

by linearizing f and h about some nominal vector solution X^*. 

This report also points out that REA's SSME model may be modified 
to form a failure detection procedure. The reader should be reminded 
that the intuitive failure detection approach outlined from equations (4.3) 
and (4.4) is based on an assumption that f and h of equations (2.1) and 
(2.2) are linear functions. Since functions f an3" h are generally non- 
linear for SSME applications, it is worthwhile to investigate the effective- 
ness of bias estimation techniques as an SSME static test failure detec- 
tion device. 
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